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Reconstruction of sparse wavelet signals from 
partial Fourier measurements 

Yang Chen, Cheng Cheng and Qiyu Sun 


Abstract —In this paper, we show that high-dimensional sparse 
wavelet signals of finite levels can be constructed from their 
partial Fourier measurements on a deterministic sampling set 
with cardinality about a multiple of signal sparsity. 

I. Introduction 

Sparse representation of signals in a dictionary has been 
used in signal processing, compression, noise reduction, source 
separation, and many more fields. Wavelet bases are well 
localized in time-frequency plane and they provide sparse 
representations of many signals and images that have transient 
structures and singularities (01, 0). In this paper, we consider 
recovering sparse wavelet signals of finite levels from their 
partial Fourier measurements. 

Let D be a dilation matrix with integer entries whose eigen¬ 
values have modulus strictly larger than one, and set M = 

| det D| > 2. Wavelet vectors = ('> • • •, V ; m,r) T 5 1 < 
m < M — 1, used in this paper are generated from a 
multiresolution analysis {Vj}j e z, a family of closed subspaces 
of L 2 := L 2 (M n ), that satisfies the following: (i) Vj C Vj+% 
for all j G Z; (ii ) V^ x = {/(D-), / G Vj} for all 
j G Z; (iii) U je %Vj = L 2 \ (iv) ^jeiMj = {0}; and (v) there 
exists a scaling vector <I> = (</>!,•• • , 0 r ) T G Vo such that 
— k), 1 < l < r, k G Z n } is a Riesz basis for Vo (0, 

0, 0, El, 0, 0, 0, 0, 0, OH). They generate a Riesz 

basis {M J / 2 ^ m (D J x—k) : 1 < m < M — 1, k G Z n } for the 
wavelet space Wj := Vj+\ © Vj, the orthogonal complement 
of Vj in Vj+ 1 , for every j G Z. Therefore any signal / in 
the scaling space Vj of level J > 0 has a unique wavelet 
decomposition, 

/ = /o + #(H- Ygj- 1, (1.1) 

where 

fo= a^(k)$(--k) e Vo (1.2) 


9j = E E b^ J (k)M^ m (D^-k) e Wj, 0 < j < J- 1. 

m= 1 k£Z n 

(1.3) 

In this paper, we consider wavelet signals / G Vj with / 0 and 
0 < j < J — 1, in the above wavelet decomposition having 
sparse representations. 
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Define Fourier transform of an integrable function / on M n 
by 

/(<£) = / /(t)e-^t. 

J R n 

Due to coherence of wavelet bases between different levels, 
the conventional optimization method does not work well to 
reconstruct a sparse wavelet signal / of finite level from its 
partial Fourier measurements /(£),£ G D, on a finite sampling 

set Q (CD, 02), H3, 02), 02), 02), d, dl). Recently, 

Prony’s method was introduced in £19), (20) for the exact 
reconstruction of one-dimensional sparse wavelet signals. 

Denote by the cardinality of a set E. We say that a 
wavelet signal / G Vj has sparsity s = (sq? • • • , sj-i) if it 
has sparsity 

ma x{#K 0 , , #K M - 1 , 0 } if j = 0 

max{#iTij,..., #K M -i,j} if j =» 1,..., J - 1, 

at level j, 0 <3 < J — 1, where Kq and K m j are supports 
of coefficient vectors (ao(k))kez™ anc ^ (bm,j( k))kez n in the 
wavelet decomposition Qi, 0 and 0 respectively. For 
the classcial one-dimensional scalar case (i.e. n = l,r = 1 
and D = 2), under the assumption that Fourier transform of 
the scaling function 0 does not vanish on (—7r,7r), 

0, £ e (—7T, 7T), (1.4) 

Zhang and Dragotti proved in C9) that a compactly supported 
sparse wavelet signal of the form 0 can be reconstructed 
from its Fourier measurements on a sampling set D of size 
about twice of its sparsity so + • • • + sj_i. In this paper, we 
extend their result to high-dimensional sparse wavelet signals 
without nonvanishing condition 0 on the scaling vector 
<f>. Particularly in Theorem |III.1[ we show that any s-sparse 
wavelet signal / of the form (|I.l|) can be reconstructed from 
its Fourier measurements on a sampling set D with cardinality 
less than 2Mr(so + ••• + sj_i), which is independent on 
dimension n. 

II. Multiresolution analysis and wavelets 

Set M = D t . Then the scaling vector = (0 1 ,..., 0 r ) T 
of a multiresolution analysis {Vj}j e % satisfies a matrix refine¬ 
ment equation, 

$(£) = G 0 (M -1 £)$(M -1 £), (III) 

where the matrix function Go of size r x r is bounded and 2tt- 
periodic. In this paper, we assume that Go has trigonometric 
polynomial entries. Hence <f> is compactly supported, and the 
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Riesz basis property for the scaling vector can be reformu¬ 
lated as that (<£(£ + 27rk))kez™ has rank r for every £ G M n . 
Therefore for any £ G M n there exist k(£, l) G Z n , 1 < / < r, 
such that 


and 

j-i 

#0 < X = 2Mr(s 0 + si + • • • + sj_i). (DI.2) 
i=o 


($(£ + 27 rk)) keAf? . has full rank r, ( 112 ) 


where 


A(0 = {k(e,Z)GZ n : (II.3) 


Let p m , 0<ra<M — 1, be representatives of Z n /MZ n , 
and write 

M—l 

z n = U (p m + MZ n ). 

m=0 

Take matrices G m , 1 < m < M — 1, with trigonometric 
polynomial entries such that 

M—l 


X Gp(£ + 27rM + 27rM- 1 p m /) =0 

(II-4) 


ra'=0 

for all 1 < m < M — 1, and 

G(£) has rank Mr for all £ G 

where 


(H.5) 


G(f) = 


/ Go(C + 27rM- 1 po) ••• Go(e + 27rM- 1 p M _ 1 ) \ 

Gi(e + 27rM- 1 p 0 ) ••• Gi(e + 27rM- 1 p M _ 1 ) 


\Gm- i($ + 27tM 1 po) ••• Gm-i (C + 27tM 1 pm-i)/ 


The following is the main theorem of this paper. 

Theorem III.l. Let D be a dilation matrix, <f> be a compactly 
supported scaling vector, 1 < m < M — 1, be wavelet 
vectors satisfying dlL4l > and B- let Ll be the set in dnu) 
with h = (hi, ..., h n ). If 1, ..., h n are linearly indepen¬ 

dent over the field of rationals, then any s-sparse wavelet 
signal of the form i p) , pt i and O can be reconstructed 
from its Fourier measurements on Ll. 


Proof: Let / be an s-sparse signal with wavelet represen¬ 
tation p) , (jOJi and Ok Set 

ao(0 = X a o(k)e“' k ' ? , (III.3) 

kez n 

and 

Vi(0= X Vi(k)e" ik€ (HI-4) 

kez n 

for 1 < m < M — 1 and 0 < j < J — 1. Then taking Fourier 
transform on both sides of the equation 0 gives 

J-l M — l 

/(0 = + X E 

j=0 m= 1 

(III.5) 

Define fi, 0 < i < J — 1, by 


In this paper, wavelet vectors T r m ,l < m < M — 1, are 
defined as follows: 

$m(0 = G m (M- 1 0$(M“ 1 0, 1 < m < M — 1. (II. 6 ) 

Then T m are compactly supported and { M :) / 2 T,,, (D y x — k) : 
1 < m < M — 1, kGZ n } forms a Riesz basis for the wavelet 
space Wj := Vj+i © V 3 for every j G Z. 

For the scaling vector <f> and wavelet vectors \k m , 1 <m< 
M — l, constructed above, one may verify that any signal in 
Vj has the unique wavelet decomposition ( [ITT] ), ( |L 2 | ) and ( [O] ). 

III. Reconstruction of sparse wavelet signals 

Take h = (hi,... ,h n ) G M n and sparsity vector s = 
(s 0 ,...,sj-i), and set |js||oo = max 0 <j<j_i sj. For 0 < 
j < J — 1 and 0 < m < M — 1 , let 

r i = {(-Si + l/2)h, (- Sj + 3/2)h, l/2)h}, 

and 

% = U 7er ,. U^“o (tt 7 + 27rM J p m 

+27rM J+1 A(7rM^- 1 7 + 2irM _:1 p m )), 


i M—l 

M) = + X E 

j =0 m=l 

(III.6) 

Then 


fj-i = f, (ni.7) 


m= 1 
M—l 

= af + E 

m=1 

M — l 

= (af(£)Go(M- 1 0+X e m,i«) 

m=l 

xG to (M -1 £))$(M -1 £) (III. 8) 

for some vectors a^(£) with trigonometric polynomial entries, 
where the last equality follows from ( |IL1| ) and ( |IL 6 | ). 

For 0 < j < J — 1, 7 G Yj and 0 < m' < M — 1, set 
% ( 7 , m') = 7 rM _J 7 + 27rp m /. 


where the set A of cardinality r is defined by ( |II.3| ). Set 

D = U jlfaj. (III.l) 


Then 


LI C {( 11s|Joo + l/2)h7r,..., (11s11oo - l/2)h7r} + 27rZ n , 


Applying ( |III.7| ) and ( |III. 8 | ) with i = J — 1, replacing £ in 
( |III. 8 | ) by 777 - 1 ( 7 , m f ) + 27rMk,k G A(M _ 1 ? 7 j_i( 7 ,m')), 
and using periodicity of aj_i and we obtain 

/(M J - 1 r; J _ 1 ( 7 ,m , )+ 27 rM J k) 

= A (J — 1,7,to / )$(M _1 t;j_i(7,to') + 27 rk) (III. 9 ) 
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for all k G A(M 1 77 . 7 - 1 ( 7 , m')), where 
A(J ~ !,7) TO ') = aj-i(7rM _J+1 7)G 0 (M _1 7?j_i(7, to')) 

M —1 

+ X] ^ )J _i(7rM _J+1 7)G m (M _1 7?j_i(7,TO')). (III.10) 

m=l 

Recall from < |II.2| > that 

f$(M- 1 r;j_ 1 (7,TO') + 27rk)) (III. 11) 

is nonsingular. Then A(J — 1 , 7 , 777 /) can be solved from the 
linear system ( |III.9|) for a ll 0 < m! < M — 1 and 7 G Tj_i. 
Recall from ( |II 4| ) and ( |II 5| ) that 


and 


G(/ittM“ J 7 ) = 


/ G 0 (M 1 777-1(7,7+)) \ 

Gi(M _1 ?7j_i( 7, m 7 )) 


(III. 12 ) 


bm,o (^^rh) — ^ v b m? o(k) (e 
k£K m 0 


V Gm-i(M 1 777-1(7, m')) ) 

is nonsingular. Thus, for every 7 G Tj_ 1 and 1 < m < M— 1, 

and Vj-^M^+by) (III. 13) 

are uniquely determined from samples of / on 12j_i C 12 by 
dnUQl ) and ( pm) . 

For 1 < m < M—1, it follows from the linear independence 
assumption of 1, hi ,..., h n on the field of rationals that 

e -z7rk-M h ^ ^ K rn ^ J _ 1? are distinct to each other. 

For 7 = nh with n G {—sj_i + 1/2,..., sj_ 1 — 1/2}, 

b m , J _ 1 (7TM- J+1 7)= £ b m , J _ 1 (k)(e-^ k - M ' J+lh ) n 

(III. 14) 

by <|III.4|». Therefore applying Prony’s method ((18), f 19j, l2H . 
G2, ED, EQ, G3, ESI ) recovers trigonometric polynomials 
b m ,j-i,l < m < M — 1, from their measurements on 
7 tM _j+1 7,7 G Tj_i. Hence b m? j_i(k), k G Z n , can be 
recovered from samples of / on 12 for all 1 < m < M — 1. 

By the above argument, 

fj-i ~ fj -2 and /j_ 2 (£), (III. 15) 

can be obtained from samples of / on 12, because 

M—l 

1-2(0 = m - E (E b m , J _ 1 (k) e - ik - M - J+1 «) 

ra=1 ke^ n 

x$ m (M- J+1 0 

by ( |III.3| ) and ( |III.8| ). Inductively we can reconstruct 

fi ~ fi-i and jU(0, (III. 16) 

from the samples of / on 12 for i = J — 2, • • • ,1. 

Taking i = 1 in ( |III.16| ) determines samples of /o on 
Next we recover the function / 0 from its Fourier measurements 
on H 0 C 12. By ( [1101 ) and pO] ), 

a 0 (n7rh) = ^ a 0 (k)(e _ * 7rk ' h )” 

keKo 


where n G {—so + 1/2, —so H- 3/2, -. -, so — 1/2}. Similar to 
( |III.13| ), we can show that 

ao(n7rh) and b m? o(/i7rh), 1 < m < M — 1, 

are uniquely determined from samples of /o on 12. Applying 
Prony’s method again recovers ao(k) for k G Ko and b m? o(k) 
for 1 < m < M — 1 and k G K m? o. Therefore /o could be 
completely recovered from its Fourier measurements on 12. 
This together with pII7| ), (IIIl 5) and pII.16[ ) completes the 
proof. ■ 

The linear independence requirement on h = (hi, . .., h n ) 
in Theorem |III.1| can be replaced by a quantitative condition 
if the sparse signal has some additional information on its 
support, c.f. na. 

Corollary III.2. Let D, <f> and 4/ m ,l < m < M — 1, be 
as in Theorem 1II.1 and let f be an s- sparse signal in CD 
satisfying 


K 0 c [a,b) n and K md C D J [a,b) r 


(III. 17) 


where 1 < rn < M — 1, 0 < j < J — 1 and a < b. Then 
f can be recovered from its Fourier measurements on 12 in 
dnuT ) with h = (h \,..., h n ) satisfying 

0 < (6 - a) (hi + h 2 + • • • + K) < 2. (III.18) 


Proof: Following the argument in Theorem III. 1 it suf¬ 
fices to prove that e _27rk ' h , k G AT 0 , are distinct, and also that 
e - zvrk-M J h ^ ^ K m j , are distinct for every 1 < m < M — 1 
and 0 < j < J — 1. The above distinctive property follows 
from pII.17| ) and pII.18| ) immediately. ■ 

From the proof of Theorem |III.1| we have the following 
result on the reconstruction of an s-sparse trigonometric poly¬ 
nomial from its samples on a set of size 2s. 

Corollary III.3. Let h = (hi,..., h n ) with 1, hi ,..., h n 
being linearly independent over the field of rationals, and 
define 

0 S = {(-s + l/2)h, (-s + 3/2)h, ...,(s- l/2)h}, 5 > 1. 
Then any n-dimensional trigonometric polynomial 


—ik 


p(o = p( k ( e 

kGZ^ 

with sparsity s, 

#{k : p(k) ^ 0} < s, 
can be reconstructed from its samples on 0 S . 

IV. Simulations 

The following algorithm for sparse wavelet signal recovery 


is proposed in the proof of Theorem III. 1 

Algorithm 1: 

1. Input sparsity vector s = (so, • • • , sj_i). 

2. Input Fourier measurements /(£),£ G 12 and set fj-i = 
/• 
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3. for j = J — 1 to 0 do 

for every 7 G Yj do 

for every m! — 0, • • • , M — 1 do 
3a) rjj ( 7 , to') = 7 rM _J 7 + 27rp TO /. 

3b) Solve the linear system 

(/.(M'Xjft.m') + 

= + 2 irk)) tsA(M _, 

to get 

■= aJ(7rM _J 7)Go(M _ 1 ^(7,TO')) 

M—1 

+ XI bX( 7 rM “ J 7)G TO (M“ 1 7 7 i (7,TO')). 

m=l 

end for 

3c) Solve the linear equation 

(aj( ttM - ^), bf ;i (ttM _j 7), • • • , 6m _ 1?5 (ttM - ^)) 
xG(kM-^7) = f4(j, 7 , ()),••• ,A(j,7,M-1)). 

end for 

3d) Recover b m?J - from b m j( 7 rM _ - 7 7 ), 7 G Ty with 
Prony’s method for every 1 < m < M — 1. 

3e) Subtract E™=i from 

/j(£) to g et fj- 1 ( 0 > £ e a 

end for 

4. Recover ao from a( 7 T 7 ), 7 G To with Prony’s method. 

5. Reconstruct the sparse wavelet signal 

/(*) = X a o( k )$(t-k) 

kez n 

J-l M —1 

+ X X X bX(k)M% m (D^t-k). 

j=0 m=l k£Z n 

Next we present simulations to demonstrate the above 
algorithm for perfect reconstruction of sparse wavelet signals 
of finite levels. Let = X[o,i)M and 02 (t) = 2y/3(t — 
l/ 2 )x[ 0 ? i)(f) be scaling functions, and let 

V’t (t) = (6 1 - l)x[o,i/ 2 )(t) + (6 1 - 5 )x[i/ 2 ,i)(*)> 

and 

Mt) = 2\/3(2f - 1/2 ) X[ o,i/2) (t) ~ 2^3(2 1 - 3/2)x[i/2,i)(*) 

be wavelet functions. Consider reconstructing the sparse signal 

/(f)=a^(2)$(f-2) + aJ(4)$(f-4) 

+ b o - 1) + bj(5)^(f - 5) 

+bf (6)^(2f - 6 ) + bf (12)^(2f - 12) (IV.l) 

from its Fourier measurements on the sampling set 

Yl = | — Y^ n7T + : n = ^ an d ^ = 0, ±1? ±2, 4 j 

in ( |III.1| ), where <f> = ( 0 i, 02 ) t , = (Vt,^) 7 ", and the 

nonzero components of ao, b 0 and bi are randomly chosen 
in [—1,1] \ (—0.1,0.1), see Figure [I] Applying the proposed 
algorithm, our numerical results support the conclusion on 




Fig. 1. Plotted on the left is the sparse wavelet signal / in ( |IV. lf while on 
the right is the magnitude of its Fourier transform and the measurements on 
ft. 


perfect recovery of sparse wavelet signals from their Fourier 
measurements on Yl. 

The proposed algorithm is tested when the Fourier measure¬ 
ments of the signal / are corrupted by random noises e, 

h(0 = f(0 + <0, 

In this case, sparsity locations obtained by Prony’s method in 
the algorithm are not necessarily integers, but it is observed 
that they are not far away from the sparsity locations of the 
signal /, when the signal-to-noise-ratio (SNR), 

SWfl = -201o Sl „ mm ^ |tKI1 

max t6 „ |/(OI 

is above 50 dB. Taking nearest integers of those locations may 
perfectly recover the sparsity positions {2,4} for the scaling 
component of level 0 , {1,5} for the wavelet component of 
level 0, and { 6 ,12} for the wavelet component of level 1. Then 
the signal / can be reconstructed by the proposed algorithm 
approximately, see Figure [2] 



Fig. 2. The difference between the original signal and the reconstructed 
signal. In this figure, the reconstruction is generated by the proposed algorithm 
with modified Prony’s method, and the noise level on Fourier measurements 
is SNR=50 dB. 


We also tested our proposed algorithm for two-dimensional 
wavelet signals with dilation D s; ^^. Presented on 

the left of Figure [3] is the amplitude of a sparse wavelet signal 

f(t i, t 2 ) = ao0(*i - 1, h) + ai0(£i - 2, t 2 - 3) 

+ 5o0(fi - 2, t 2 - 1) + 5i0(£i - 3, t 2 - 5), (IV.2) 

where ao,ai,&o>fri £ [—1,1] \ (—0.1,0.1) are se¬ 
lected randomly, the scaling function is 4>(ti,t 2 ) = 

X[o,i)(^i)X[o,i)(^ 2 )> and the wavelet function is ^(ti,t 2 ) = 
X[o,i)(^i)(X[o,i/ 2 )(^i) - X[i/ 2 ,i)(* 2 ))- Oar simulations show 
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that the signal / in ( |IV.2| ) can be reconstructed from its Fourier 
measurements on 


n 


{{^§ n+2k ^ n+2l h 


= d=l, ±3 and k,l = 0,1 J>, 
(IV. 3) 

which is plotted on the right of Figure [3] 




Fi g. 3. Fourier amplitudes of the signal / in pV.2) and the sampling set £2 
in ( |IV.3} for sparse recovery. 


V. Conclusion 

In this paper, we show that sparse wavelet signals of finite 
level can be reconstructed from their Fourier measurements on 
a deterministic sampling set, whose cardinality is independent 
on signal dimension and almost proportional to signal sparsity. 
A difficult problem on this aspect is exact reconstruction 
of signals having sparse wavelet-like (e.g. wavelet packet, 
framelet, curvelet, and shearlet) representations from their 
partial Fourier information ([27), EEL [29)> 1 301 . ED ). 
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